Temperature dependence of spin-polarized transport in ferromagnet / unconventional 

superconductor junctions 
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Tunneling conductance in ferromagnet / unconventional superconductor junctions is studied tlie- 
oretically as a function of temperatures and spin-polarization in ferromagnets. In d-wave supercon- 
ductor junctions, a zero-energy Andreev bound state drastically affects the temperature dependence 
of the zero-bias conductance (ZBC). In p-wave superconductor junctions, numerical results show var- 
ious temperature dependence of the ZBC depending on the direction of the magnetic moment in 
ferromagnets and the pairing symmetry in superconductors such as p^, Py and Pa, -I- ipy-w&ve sym- 
metries. The last one is a candidate for the pairing symmetry of Sr2Ru04. From these characteristic 
features in the conductance, we may obtain the information about the degree of spin-polarization 
in ferromagnets and the direction of the d-vector in spin-triplet superconductors. 



I. INTRODUCTION 

In recent years, transport properties in unconventional 
superconductor junctions have been studied both the- 
oretically and experimentally. In these junctions, a 
zero energy state (ZES)^'-'^ formed at the jimction in- 
terface plays an important role in the tunneling spec- 
troscopy. It is now well known the ZES is responsible 
for zero-bias conductance peak (ZBCP) in the high-Tc 
superconductor junctions^iS*SiLSiS and related phcnom- 
3^^io,ii,i2,i3,i4,i5,i6,i7,i8,i9,2Q,2i_ ^he theoretical studies 

clearly relate the formation of the ZES to the ZBCP 
in tunneling spectroscopji2Si2i2ii2^. Since the formation 
of the ZES is a general phenomenon in unconventional 
superconductor junctions, the ZBCP is also expected 
in spin-triplet superconductor junctions^SiSLSSiSSi^SiSiiS. 
Actually, the ZBCP has been observed in junctions 
of SraRuO^^^ and UBci^S. The ZBCP has also 
been theoretically predicted for organic superconductors 
(TMTSF)2X very recently ?'Pi?7i??. 

From a view of future device application, transport 
properties in hybrid structures consist of ferromagnets 
and superconductors have attached much attention. It 
was pointed out in ferromagnet / insulator / spin-singlet 
unconventional superconductor (F/I/S) junctions that 
the amplitude of the ZBCP decreases with increasing 
the magnitude of the exchange potential in ferromagnets. 
This is because the exchange potential breaks the time- 
reversal symmetry and suppresses the retro-reflectivity 
of the Andreev reflection 39,40,4i,42,43,44,45,46,49,5Q^ Thus 
the ZBCP is sensitive to the degree of spin-polarization 
in ferromagnets. Since the tunneling conductance is 
independent of the magnitude of the insulating bar- 
rier—, it is possible to estimate the spin-polarization 
in ferromagnets through the temperature dependence of 
the ZBCP. An experimental test would be carried out 
in Lao.7Sro.3Mn03/YBa2Cu307_i: junctions in near fu- 
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When spin-triplet superconductors are attached to fer- 



romagnet 



,43.44 



the ZBCP depends not only on the spin- 



polarization also on other parameters such as relative an- 
gles between the d-vector in triplet superconductors and 
the magnetic moment in ferromagnets2i4S4Si^. Thus it 
may be possible to know details of the pair potential by 
comparing the characteristic feature of the ZBCP in the- 
oretical calculations and those in experiments. For this 
purpose, it is necessary to know efl'ects of another ingre- 
dients such as temperatures and the profile of the pair 
potential near the junction interface on the ZBCP. It is 
known that the amplitude of pair potential is drastically 
suppressed at a surface or a interface of superconduc- 
tors in the presence of the ZES2.28,56,57j58^_ Although 
there are several studies on tunneling phenomena in fer- 
romagnet / unconventional superconductor junctions so 
fari^^i^, such issues have never been addressed yet. 

In this paper, we calculate the tunneling conductance 
in ferromagnet / unconventional superconductor junc- 
tions as a function of temperatures and degrees of spin- 
polarization in ferromagnets, where the spatial depen- 
dence of the pair potential is determined self-consistently 
based on the quasiclassical Green's function theory. We 
choose d-wave, and Px + ipy-wave symmetries for the pair 
potentials which are candidates for pairing symmetries of 
high-T(7 cuprates and Sr2Ru04, respectively. For com- 
parison, we also study the conductance in p^- and Py- 
wave superconductor junctions. From the calculated re- 
sults, we reach the following conclusions. 

(1) In d-wave junctions with (110) orientation and p^- 
wave junctions, an incident quasiparticle from a ferro- 
magnet always feels the ZES irrespective of the incident 
angles. The zero-bias tunneling conductance (ZBC) at 
the zero temperature is insensitive to the barrier poten- 
tial at the interface. This result indicates a possibility to 
estimate the magnitude of the spin-polarization of ferro- 
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magnets at sufSciently low temperatures (T) in experi- 
ments. 

(2) In d-wave junctions with (110) orientation and Px- 
wave junctions, the ZBC monotonically decreases with 
increasing temperatures for small magnitudes of spin- 
polarization. On the other hand for large magnitudes 
of polarization, the ZBC becomes an increasing function 
of temperatures. While for d-wave junctions with (100) 
orientation and Py-wave junctions, where the ZES does 
not appear, the ZBC is an increasing function of T in- 
dependent of the spin-polarization. For p^ + ipy-wave 
junctions, the ZBC first decreases with increasing T then 
increases. 

(3) For p-wave junctions, the ZBC has various tempera- 
ture dependence depending on the direction of the mag- 
netic moment in ferromagnets. This unique property is 
peculiar to spin-triplet superconductors. 

(4) Throughout this paper, we calculate the ZBC in two 
ways; i) the spatial dependence of the pair potential is 
assumed to be the step function (non-SCF calculation), 
ii) the spatial depletion of the pair potentials is deter- 
mined self-consistently (SCF-calculation). By compar- 
ing the conductance in the two ways, we found that the 
results in the non-SCF calculation are qualitatively the 
same with those in the SCF-calculation. 

The organization of this paper is as follows. In Sec. 2, 
we formulate the tunneling conductance with arbitrary 
angle between the magnetization axis of the ferromagnet 
and c-axis of the superconductor. We show the tunneling 
conductance depends on the direction of the magnetic 
moments only when superconductors have spin-triplet 
pairing. In Sec. 3, we calculate the polarization and tem- 
perature dependence of ZBC for both d-wave and p-wave 
junctions. In Sec. 4, we summarize this paper. 



H{r)S{r - r' 
At(r,r') 



A(r,r') 

-H*{r)6{r-r') 



where E is the energy of a quasiparticle measured from 
Ef, H{r) = /ioi-U(r)-^(r), ho = -^V^ + Vix)- Ep, 
U(r) = UQ{—x)n, 1 and <t are the 2x2 identity matrix 
and the Pauli matrix, respectively. Here n points the 
direction of the magnetic moment in ferromagnets and 
<d{x) is the Heaviside step function. The indices t and J, 
denote the spin degree of freedom of a quasiparticle in su- 
perconductors. The magnetization axis in ferromagnets 
is represented in a polar coordinate {9mt4>m) ss shown 
in Fig. 1. We assume that the quantization axis of spin 
in the triplet superconductors is in the c-axis which is 
parallel to z direction. At first, we assume that the pair 
potential is a constant independent of x, 



A{es,x)^^{esmx), 



(2) 



where /S.{9s) does not have the spatial dependence, kx — 
kp cos 9s and ky = kp sin6s are the wavenumber in su- 
perconductors with kp being the Fermi wave number. 
The pair potential A{Os) is given by 



Axt(^s) Au(^?s) 



(3) 



The Hamiltonian in Eq. (|2Jl is written in the coordinate 
of the spin space in the superconductor. It is compre- 
hensive to rewrite the Hamiltonian in the coordinate of 
spin space in the ferromagnet since this notation is useful 
to consider the scattering processes. The Hamiltonian in 
the coordinate of spin space in the ferromagnet is ob- 



II. FORMULATION 



Let us consider a two-dimensional F/I/S junction in 
the clean limit as shown in Fig. ^ We assume a flat 
interface at x = 0. The insulator is described by the 
delta-function V{x) = H6{x), where H represents the 
strength of the barrier potential. We also assume that 
the Fermi energy Ep and the effective mass m in the 
ferromagnet are equal to those in the superconductor. 
The Stoner model is applied to describe ferromagnets, 
where the exchange potential U characterizes the ferro- 
magnetism. The wave numbers in the ferromagnet for 
the majority (t) and the minority (|) spins are denoted 

The wave functions ^'(r) 



hyk-^a)-^p{Ep + i-)U). 

are obtained by solving the Bogoliubov-de Gennes (BdG) 
equation under the quasiclassical approximationSiffl 



E^{r) = / dr'J?(r,r')*(r') , *(r) = 




(1) 





U j 








^ 



ferromagnet 



superconductor 



FIG. 1: Schematic illustration of a ferromagnet / supercon- 
ductor junction. The direction of the magnetization axis is 
denoted by a polar coordinate {Om,4im)- 



3 



tained by using the following unitary transformation: 



U = 



jj 

IJ- 



U = 
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72 
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(4) 
(5) 

(6) 



where U is the operator which diagonalizes the H{r). 
The effective pair potential in the coordinate of spin 
space in ferromagnet is rewritten as 



(7) 



Here, we consider the four types of pair potentials, d- 
wave, Px-'wave, py-wave and Px + ipy-wave symmetries 
in superconductors. In d-wave case, the pair potential is 
described as 

Aui9s)^-A^^i9s)^Aof{es), (8) 
A„(05)=Aa(^^s) = O, (9) 
f{9s)=cos[2ies-a)], (10) 

where a is the angle between a-axis of the high-Tj, super- 
conductors and the interface normal. The effective pair 
potential in the coordinate of spin space in ferromagnet 
is given by 



(11) 

(12) 
(13) 




Aof{0s) 
^Aofies) 



=A{9s). 



As shown in Eq. H13I) , the expression of the pair potential 
remains unchanged under the transformation in Eq. |0J. 
Therefore transport properties are expected to be inde- 
pendent of the direction of the magnetic moment. This 
conclusion can be applied to any spin-singlet supercon- 
ductors. On the other hand in spin-triplet superconduc- 
tors, the pair potentials are given by 



AofiOs), 
0, 



(14) 
(15) 



where the direction of the d-vector is parallel to the c-axis 
and 

{cos 9s for Px — symmetry, 
sin 6*5 for - symmetry, (16) 

qiBs fQj. _ symmetry. 

Because the spin degree of freedom of Cooper pairs is 
active in triplet superconductors, the pair potential after 
the transformation in Eq. Q depends on the direction 
of the magnetic moment 



sin 9 m cos 9 m 
cos 9]^J — sin 9 m 



f{9s)Ao. (17) 



There are four reflection processes when an electron with 
the majority spin is incident from ferromagnets: 

i) Andreev reflection to majority spin (a^^) 

ii) Andreev reflection to minority spin (a^i) 

iii) normal reflection to majority spin (5||) and 

iv) normal reflection to minority spin (6|j). 

Similar reflection processes are also possible, when an 
electron with the minority spin is incident from ferromag- 
nets. The Andreev and the normal reflection coefficients 
are denoted by aggi and bge' , respectively. In these coef- 
ficients, a quasiparticle is reflected from the spin-channel 
s into the spin-channel s' . 

The wave function in ferromagnets for majority spin 
injection is represented by 




1 





Vl 



/O 




(18) 



where kpj < ks < kpt and kf 



s^p^ Ks <- Kp^ ana ks ~ y • The wave 
function for minority spin injection is written in the sim- 
ilar way. The coefficients agg' and bggi are determined by 
solving the BdG equation with the quasiclassical approx- 
imation under appropriate boundary conditions. 

The tunneling conductance axieV) for finite temper- 
ature is given bjiSLSSiS 



2e2 



G = 



1 



16/cbT 

X (c^sf(^s) + (^sli^s)) sech' 



dE I d9s cos 9s 

■it/2 

2 fE~eV 



2kBT 



(19) 



(20) 



o-sf =1 + I«TtI^ - 1^ 



I 



^1 



ni 



|2_!ZI 15-12 



^Q{9c-\9s\), 



(21) 



'^51 = 1 i + ?^NtP + Nil'- - 



xQ{ec-\9s\). 

with Zg^ = Z/ cos 05, Z = 
v/l± A:/cos2 6*5. Here X = 



(22) 



2mH/h^kp and ?7|(j) — 

U I Ep is defined as the 
spin-polarization parameter. The quantity o'sf(j) is the 
tunneling conductance for an incident electron with the 
majority (minority) spin. For \9s\ > 9c = cos~^ VA", the 
refiected wave becomes an evanescent wave and does not 
contribute to the tunneling conductance. As shown in 
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above equations, the tunneling conductance depends on 
6m only when superconductors have spin-triplet Cooper 
pairs. The tunneling conductance can be summarized in 
simple equations in following several cases. When 9m is 
or TT in spin-triplet superconductors, the conductance 
is described by 



A=[i-|r+r_|^(i-a^j)-fa^j|r+|1 

X Q{ec-\Os\)/LDu 

B^i- Q{ec - \esm- I r+r_ \'']/Ldi, 
c=[i-|r+r_|2(i-a^^) + a^^|r+p] 

X <d{ec-\0s\)lLD2, 



L 



Dl 



Lni — 



i-r+r_ji 



1 - cr 



X exp [i{Lpi - if-^)] 

1 - r+r 



with 



X exp - ifi)] 1^ , 
exp(iv5j) 



1 



iZa 



exp(-iv5|] 



^l-cr^j(l + ?7j -iZos)' 



(23) 
(24) 



(25) 
(26) 

(27) 
(28) 
(29) 

(30) 
(31) 



Above equations can be applied to the spin-singlet su- 
perconductors. When 9m — 7''/2 in spin-triplet super- 
conductors, the conductance for each spin is given by 



i- |r+r_p(i-a^^) + a^^|r+|^ 
|i-r+r_(i-a^^)P 
|r+r_|2(i-a^j) + a^i|r+|2 



|i-r+r_(i-a^i) 
xe{dc-\es\). 



In above equations, we define 



At] 



'Nl 



(l+r;j)2 



ei0c-\es\). 



(32) 

(33) 

(34) 
(35) 
(36) 



In this paper, the dependence of the pair potential on 
temperatures is described by the BCS's gap equation. 
The spatial dependence of the pair potential can be de- 
scribed by A{&S,x) with A(6's,a;) = A{0s,x)e{x). In 
order to determine the spatial dependence of A{9s,x), 



we apply the quasiclassical Green's function theory de- 
veloped by Hara, Nagai, et. al. 2,58^64,65 ^ ^^le following, 
we briefly explain the method in the case of spin-singlet 
superconductors. An extension to spin-triplet supercon- 
ductors is straightforward. The spatial dependence of 
A (0s, a;) is calculated by the diagonal elements of the 
matrix Green function gaa{x) which are represented by 



l + D+{x)F+ix) 



2iF+(x) 



n,,(fj^T\-n ^-D+ix)F+{x) l-D+(x)F+{x) forj) 

9++\.^S,X) -I I 2iD+(x) l+D^(x)F+(x) I ' l-J'i 

\ l-D^{x)F+{x) l-D+{x)F+(x) J 

/ l+D-{x)F^{x) 2iF^ jx) 



-l+D_(x)F_(x) -l+D-{x)F-{x} 
2iD-(x) l+D-(x)F-(x) 
-l+r>_(x)F_(x) -l+D^{x)F-{x) 



(38) 



where an index a = ± specifies the direction of the mo- 
mentum in the x direction. In these Green functions. 
Da (x) and Fa (x) obey the following equations 

^VFx\Da{x) = a 

X [2lu„Mx) + A{es,x)Dl{x) - A*{es,x)] , (39) 

^VFx\Fa{x) = a 

X [-2L0mFa{x) + A*{es.x)Fl{x) - A{es,x)] , (40) 
MOs^x) 



A{es,x) = 



-Aies,x) 







(41) 



The boundary conditions at the interface are given 
by^MI 

(fi-. - I + iZ)(m - I - iZ) 
F+W =T^ 7Tp^ -(d_(0)-\ (42) 



(t]-^ + 1 + iZ){r]i + 1 - iZ) 
{r,-^-l + iZ){r,i-l-iZ) 
{f]-^ + l + iZ){r]i + l-iZ) 



D+{0). (43) 



We first solve D±{x) and F±{x) in Eqs. ^ and 
then calculate g±^±{0, x) in Eqs. IpTfjl and for a given 
A{es,x). By using g±,±{9,x), A{es,x) is given by the 
following equations 



.7r/2 

A(%,x)=V/ d0s'V{es,e's)ga,a{^'s,x), (44) 

v{e,e') ^gofi0)f{e'), 

2-KkBT 



9o 



'HT/Tc) + Eo 



(45) 
(46) 



<m<m„ m+1/2 



where we introduce the cut-off ma to regularize go and 
f{6) is given in Eqs. 1)1 0|l and 1)16(1 . The iteration is car- 
ried out until the sufhcient convergence is obtained. In 
this way, we obtain A{9, x), i.e., A{9, x) self-consistently. 

Under the pair potential in the self-consistent calcula- 
tion, we obtain r± (x) by solving 

ih\vFx\t+{x) 

=a [2Er+{x) ~ A{0s,x)tl{x) - A*{0s,x)] , (47) 
ih\vFx\i'-ix) 

=a [2Er_{x) - A*ies,x)r^_{x) - Aies,x)] . (48) 



5 



By substituting r± into Eqs. (|23|I - H34() . we can calculate 
the tunneling conductance. In what follows, we assume 
that the transition temperature of ferromagnets is much 
larger than Tc which is the transition temperature of 
superconductors. In such situation, we can neglect the 
temperature dependence of X. 



III. RESULTS 



polarization dependence of zero-bias 
conductance 



In this subsection, we show calculated results of the 
ZBC at the zero temperature as a function of the spin- 
polarization in ferromagnets {X). The dimensionless 
ZBC (r) is given by 



r = 



(7TiO)h 

2e2 ■ 



(49) 



At first we show the conductance obtained in the step- 
function model, where the spatial dependence of the pair 
potential is not determined self-consistently (non-SCF 
calculation). The X-dependence of F at the zero tem- 
perature for d-wave junctions is plotted in Fig. 2. The 
magnitude of F is always a decreasing function of X. 
For a — (Fig. 2(a)), F decreases with increasing of Z. 
On the other hand, for a = tt/A (Fig. 2(b)), F is com- 
pletely independent of Z because of the perfect Andreev 
reflection due to the zero-energy resonance state at the 
interface. 



e/jr 



FIG. 2: X dependence of the zero-bias conductance F in 
non-SCF calculation for (a.) a — and (b) a = 7r/4 in d-wave 
junctions at zero temperature, a; Z — 0, h: Z = 1 and c: 
Z ^5. 



Secondly, we show the polarization dependence of F 
in triplet p^, Py and px + ipy-wave superconductor junc- 
tions as shown in Fig. 3, where the pair potential are 
given in Eqs. (|14|l . 1)15(1 and ((16|l . The conductance de- 
pends on 9m for all pairing symmetries. The spin de- 
gree of freedom remains in spin-triplet superconductors. 
As a consequence, the conductance depends on the rela- 
tive angle between the magnetic moment in ferromagnets 
and d-vector in superconductor. This is the characteris- 
tic feature of p-wave junctions. For 6m = 0, as shown 



in Figs. 3 (a), (d), and (g), F approaches to zero in the 
limit oi X ^ 1 independent of Z. In these cases, the 
diagonal elements in Eq. Q disappear and a quasipar- 
ticle suffers the spin- flip in the Andreev reflection, (i.e., 
a| I = aj-^ = 0). For 16*51 > 6c, the Andreev reflection to 
I spin becomes the evanescent wave. At the same time, 
an incident wave with J, spin vanishes. Thus F vanishes in 
the limit of X — 1, where ferromagnets are referred to as 
half-metals. On the other hand, for sin6'M 7^ 0, F takes 
finite values even in X ^ 1 as shown in Figs. 3 (b) (c), 
(e) (f) (h), and (i) because the spin-conserved Andreev 
reflection is still possible in these junctions. The results 
for dM — Ti'/Z are shown in Figs. 3 (c) (f) and (i). In this 
case, the off-diagonal elements in Eq. (I17|l become zero. 
Thus spin of a quasiparticle is conserved in the Andreev 
reflection, (i.e., j = a| j = 0). The Andreev reflec- 
tion of an electron with "f spin survives irrespective of 
6s, whereas that of an electron with | spin vanishes for 
l^sl > ^c- We note in p^^-wave junctions that F does not 
depend on Z as well as d-wave junctions with a = n/A. 
This is because the ZES's are formed at the interface. 




b 


(d) 


c 




a 


(e) 


b 




c 


a 


(f) 


b 




c 




10 10 

polarization paiameter X 

FIG. 3: X dependence of r in non-SCF calculation, (a.) Om — 
0, (b) 9 At = 7r/4 and (c) 9m ~ 7!"/2 for pa;-wave junctions, 
(d) 6m — 0, (e) 9m — 7r/4 and (f) 9m — 7r/2 for p^-wave 
junctions, (g) 9m = 0, (h) 9m ~ 7r/4 and (i) 9m ~ 7i"/2 for 
Px + ipy-wsve junctions, a: Z = 0, h: Z = 1 and c: Z = 5. 



Thirdly we show the tunneling conductance under the 
pair potential whose spatial dependence is determined 
self-consistently ( SCF calculation). The results for d, 
Px, Py, and Px + ipy-wsLve junctions are shown in Figs. 4 
and 5. The conductance in SCF calculation in Figs. 4 and 
5 should be compared with corresponding results in non- 
SCF calculation in Figs. 2 and 3, respectively. We do 
not find any remarkable differences between the results 
in SCF calculation and those in non-SCF calculation as 
shown in these figures. 

Finally, X dependence of F is plotted for various 
temepratures. As seen from Fig. 2(b), using tunnel- 
ing through ZES, we can determine the magnitude of X 
thtough the value of F. This is a unique propetry for d- 
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accesible in the experiments. 



polarization parameter 



FIG. 4: X dependence of F in SCF calculation for d-wave 
junctions, (a) a = and (b) a = tt/A. a: Z = 0, h: Z = 1 
and c: Z = 5. 




10 10 1 

polarization parameter X 

FIG. 5: X dependence of F in SCF calculation, (a) 9m = 
and (b) 8m = 7r/2 with p^^-wave junctions, (c) 6m = and 
(d) 9m = 7r/2 with p^-wave junctions, (e) 9m = and (f) 
9m = 7r/2 with px + ipy-wave junctions, a: Z = 0, h: Z = 1 
and c: Z — 5. 



wave superconductor with a = 7r/4, or p-wave junctions 
where all quasiparticles feel ABS independent of their 
directions of motions. Since F is plotted at zero tem- 
perature in Fig. 2(b), it is actually important how this 
property holds even in finite temepratures. As shown in 
Fig. 6, when the magnitude of the temperature T is suf- 
ficiently smaller than Tc, F is almost insensitive to the 
magnitude of Z and we can estimate the magnitude of X 
through F. In the actual experiments, high Tc cuprates, 
e.g., YBaCuO, BiSrCaCuO, with (110) oriented interface 
is a promising candidate. In such a case, the actual value 
of O.OOlTc becomes O.IK and this temperature is fully 




polarization parameter X 



FIG. 6: X dependence of F for various temperatures for 
d-wave junctions for a = 7r/4 in non SCF calculation, (a) 
T/Tc = 0.001, (b) T/Tc = 0.005 and (c) T/Tc = 0.01. a: 
Z ^0,h: Z = 1 and c: Z ^ 5. 



B. Temperature dependence of zero-bias 
conductance 



In this subsection, we discuss the temperature depen- 
dence of F. In the first part, we show the conductance in 
the non-SCF calculation. Then the results are compared 
with those in the SCF calculation in the second part. At 
first let us focus on d-wave junctions for a = as shown in 
Figs. 7 (a), (b), (c), where the conductance is plotted as 
a function of temperatures for several magnitudes of the 
exchange potential X. We note in these junctions that 
the ZES is not formed at the interface. For Z — (see 
Fig. 7(a)), the exchange potential in ferromagnets sig- 
nificantly affects the temperature-dependence of F. For 
large X, F increases with the increase of T as shown in 
the curve d in Fig. 7(a). The Andreev reflection (two- 
electron process) is suppressed by the large exchange po- 
tential and the current is mainly carried by single elec- 
tron process. While for small X, F decreases with the 
increase of T as shown in curve a in Fig. 7(a). This is 
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because the current at the zero- voltage is mainly carried 
by two-electron process through the Andreev reflection 
and the amplitude of the Andreev reflection is suppressed 
for T Tc- For Z = 5 (Fig. 7(c)), F becomes small 
around T ~ because the insulating barrier suppresses 
the Andreev reflection. The results show that F increases 
monotonically with increasing temperatures independent 
of X since the current is mainly carried by single electron 
process. For Z — 1 (Fig. 7(b)), excepting for large X, 
the magnitude of F has a non-monotonic temperature de- 
pendence, since the amplitude of single-electron process 
(Andreev reflection) is enhanced (suppressed) with the 
increase of T. 
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FIG. 7: Temperature dependence of F in non-SCF calculation 
for d- wave junctions, (a) Z — 0, (b) Z = 1 and (c) Z = 5 with 
a = 0. (d) Z = 0, (e) Z = 1 and (f) Z = 5 with a = 7r/4. a: 
X = 0, b: X = 0.7, c: X = 0.9 and d: X = 0.999. 



Secondly we show the temperature dependence of F in 
d-wave junctions with a = 7r/4 in Figs. 7(d), (e), and (f). 
For Z = 0, the line shape of the all curves in Fig. 7(d) 
are qualitatively similar to those with a = shown in 
Fig. 7(a) since ZES is not formed at the interface. We 
note in Figs. 7 (d), (e) and (f) that F at the zero tem- 
perature is independent of Z . The Andreev reflection is 
perfect in the limit of Z = 0. For finite Z, the ZES is 
formed at the interface which also leads to the perfect 
Andreev reflection at T = 0. The characteristic behavior 
of the resonant tunneling can be seen in the conductance 
for large Z . The amplitude of F is proportional to the in- 
verse of T for intermediate temperatures (curve a in Fig. 
7 (f)). Since the retro-reflectivity of Andreev reflection 
is broken by the exchange potentials, the degree of the 
resonance at the interface is weakened. As a results, the 
temperature dependence deviates from the inverse of T 



in curbs b, c and d in Fig. 7 (f). For sufficiently large 
magnitudes oi X,T becomes an increasing function of T 
for T > 0.5Tc as shown in Fig. 7(f). In the case of the 
d-wave with a — 7r/4, we can estimate the magnitude of 
X from the temperature dependence of F. 

Thirdly we show temperature dependence of F in p^- 
wave junctions with Om = in Figs. 8(a), (b), (c) and 
those with tt/2 in Figs. 9(a), (b), (c). As shown in 
Figs. 8(a), (b) (c), the temperature dependences of F are 
very similar to those of d-wave junctions with a = 7r/4. 
When 6m = 0, F for small X is a decreasing function 
of T, whereas F for large X is an increasing function 
of T. In the case of 9m — 7r/2, however, F becomes a 
monotonic decreasing function of T independent of Z as 
shown in Figs. 9(a), (b), (c). The spin of a quasiparticle 
is always conserved in the Andreev reflection in this case. 
Therefore the suppression of the conductance due to the 
breakdown of the retro-reflectivity becomes weaker than 
that in the case of 6m = 0. 

Next we show the conductance in py-wave junctions 
with 6m = in Figs. 8(d), (e), (f) and those with n/2 
in Figs. 9(d), (e), (f). In dM = 0, the temperature 
dependence of F are similar to those of d-wave junctions 
with a = shown in Figs. 7(a),(b),(c). In these junctions, 
no ZES is expected at the interface. In 6i\i = tt/2 as 
shown in Figs. 9(d), (e), (f), F for large X are larger 
than those with dM = in low temperatures. In the 
case of = 7r/2, the spin of a quasiparticle is conserved 
in the Andreev reflection, therefore, the suppression of 
conductance due to the breakdown of retro-reflectivity 
becomes weaker than that in the case of 6m = 0. 

Finally we show the conductance in -I- ipj^-wave junc- 
tions with 6m = in Figs. 8(g), (h), (i) and those with 
7r/2inFigs. 9(g), (h), (i). Inp3;-|-tpy-wave junctions with 
6m = 0, the temperature dependence of F can be under- 
stood by the combination of the results in p^-'waye and 
those in py-wave junctions because the ZES is only ex- 
pected for a quasiparticle incident perpendicular to the 
interface?-. As seen from Fig. 8 (g) for Z~Q, there is 
no clear difference between T in px + ipy-wave junctions 
with 6m = and corresponding results in p^ or py-wave 
junctions shown in Figs. 8(a) and (d). For a finite bar- 
rier potential at Z — 1, the line shape of the all curves 
in Fig. 8(h) is rather similar to corresponding results in 
the px-wave junctions than those in the Py-wave junc- 
tions. However, F are smaller than those in the p^-wave 
junctions. For Z — 5, T for small X are enhanced at 
low temperatures because of the ZES [see curve a in Fig. 
8(i) ]. On the other hand, for large X, F is an increasing 
function of T [see curve d in Fig. 8(i) ]. For 6m — 7r/2, F 
becomes a decreasing function of T for all cases as shown 
in Figs. 9(g), (h), (i). These features are similar to those 
in the Pa;-wave junctions. 

It is important to check how the above results are mod- 
ified in the presence of the spatial dependence in the pair 
potential near the interface. We show the conductance in 
SCF calculation in the second part of this subsection. In 
what follows, we consider two cases oi X, X = (curve 
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FIG. 8: Temperature dependence of F in non-SCF calculation 
for Qm = 0. (a) Z = 0, (b) Z = \ and (c) Z = 5 in p^j-wave 
junctions, (d) Z — 0, (e) Z = 1 and (f) Z = 5 in py-wave 
junctions, (g) Z = 0, (h) Z = 1 and (i) Z = 5 in px + ipy- 
wave junctions, a: X = 0, h: X = 0.7, c: X = 0.9 and d: 
X = 0.999. 
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FIG. 9: Temperature dependence of F in non-SCF calculation 
for Om = tt/2. (a) Z = 0, (b) Z = 1 and (c) Z = 5 in p^j-wave 
junctions, (d) Z = 0, (o) Z = 1 and (f) Z = 5 in py-wave 
junctions, (g) Z = 0, (h) Z = 1 and (i) Z = 5 in px + ipy- 
wavo junctions, a: X = 0, b: X = 0.7, c: X = 0.9 and d: 
X = 0.999. 



a) and X = 0.9 (curve b). The corresponding results in 
non-SCF calculation are curves a (X=0) and c (X=0.9) 
from Figs. 7 to 9. 

In Fig. 10, we show the conductance in SCF calculation 
in the d-wave junctions with a = and 7r/4. It is found 
that the temperature dependences of F with a = in the 
SCF calculation arc very similar to those in the non-SCF 
calculation when we compare the results in Figs. 7(a), 
(b), (c) with those in Figs. 10(a), (b), (c). The same ten- 
dency can be seen between the results the d-wave junc- 
tions with a = 7r/4 in Figs. 7 (d), (e), (f) and those in 



Figs. 10(d), (e), (f). When the ZES arc formed at the in- 
terface, the profile of the pair potential significantly devi- 
ates from the step-function. The characteristic behavior 
of the conductance in the SCF, however, is qualitatively 
the same with those in the no-SCF. From the calculated 
results, wc conclude that the conductance is insensitive 
to the profile of the pair potential. This is because the 
resonant tunneling through the ZES dominates the con- 
ductance. We note that the ZES is a consequence of the 
sign-change of the pair potential. 



a=0 




1 1 
normalized temperature T/T ^ 

FIG. 10; Temperature dependence of F in SCF calculation 
for d-wave junctions, (a) Z = 0, (b) Z = 1 and (c) Z = 5 
with a = 0. (d) Z = 0, (e) Z = 1 and (f) Z = 5 for a = 7r/4. 
a: X = and b: X = 0.9. 



In the Px-wave junctions, line shapes of F for 0m = 
shown in Figs. 11(a), (b), (c) are very similar to those in 
the d-wave junctions with a = 7r/4. For 6m = 7r/2, as 
shown in Figs. 12(a), (b), (d), the magnitudes of F are 
slightly larger than those in Figs. 11(a), (b), (c). These 
features arc almost similar to those found in non-SCF 
calculation in Figs. 8(a), (b), (c) and Figs. 9(a), (b), (c). 

As well as in the Px-wave junctions, the characteris- 
tic behavior of F in SCF results in pj,-wave junctions 
shown in Figs. 11(d), (e), (f) and Figs. 12(d), (e), (f) are 
almost the same with those obtained in non-SCF calcu- 
lation shown in Figs. 8(d), (e), (f) and Figs. 9(d), (e), 
(f). 

In the case of px + ipy-wavc symmetry, the line-shapes 
of F for (g), (h) and (i) in Figs. 11 and 12 are simi- 
lar to those in non-SCF results shown in (g) and (h) in 
Figs. 8 and 9. However, the temperature dependencies of 
r based on the SCF calculation deviate from those in the 
non-SCF one for large Z [Fig. ll(i) and Fig. 12(i)]. In 
the SCF calculation, F first decreases with the increase of 
T then increases. The decreasing part is similar to that 



9 




normalized temperature T/T^- 



FIG. 11: Temperature dependence of F in SCF calculation 
for 6m = 0. (a) Z = Q, (b) Z = 1 and (c) ^ = 5 in pa;-wave 
junctions, (d) Z — 0, (e) Z = 1 and (f) = 5 in py-wave 
junctions, (g) Z = 0, (h) Z — 1 and (i) Z = 5 in paj + ipy-wave 
junctions, a: X = and h: X = 0.9. 



p, p. 




10 10 1 



normalized temperature T/T ^ 

FIG. 12: Temperature dependence of F in non-SCF calcula- 
tion for Om = 7i-/2. (a) Z = 0, (h) Z ^ 1 and (c) Z = 5 in 
Pii-wave junctions, (d) Z — 0, (e) Z = 1 and (f) = 5 in 
Py-wave junctions, (g) Z — 0, (h) Z = 1 and (i) Z — 5 in 
Pa: + ipy-wave junctions, a: X = and b: X = 0.9. 



in the p^-'waye case and is explained by the ZES. The in- 
creasing part is similar to that in the py-wave junctions. 
Since only a quasiparticle injected perpendicular to the 
interface contributes to the ZES, the effects of the spatial 
dependence of the pair potentials on the conductance are 
not be negligible^i. 



IV. SUMMARY 

In this paper, we have calculated the polarization and 
temperature dependence of the zero-bias conductance 



(ZBC) in F/I/S junctions, where we have chosen the 
symmetry of the pair potential as d-wave for high Tc 
cuprates and Px + ipy-w&Ye for Sr2Ru04. As a reference, 
we have also studied the conductance in px and py-wave 
junctions. We have established a formalism of the ZBC 
which is available for the arbitrary 6^ which is the an- 
gle between the magnetization axis in ferromagnet and 
c-axis of superconductors. The 9m dependence of the 
tunneling conductance only appears in the spin-triplet 
superconductor junctions. On the basis of the numerical 
results, we reach the following conclusions. 

(1) When injected quasiparticles from ferromagnets al- 
ways feel the zero-energy resonance state, e.g., d-wave 
junction with (110) orientation and pa;-wave junction, 
the zero-bias tunneling conductance (ZBC) at the zero 
temperature is insensitive to the barrier potential at the 
interface. This property is useful for the determination 
of the degree of the spin-polarization in ferromagnets at 
sufficiently low temperatures. One of the promising can- 
didate is LaSrMnO/YBaCuO(BiSrCaCuO) with well ori- 
ented (110) oriented interface. Within our theory, below 
O.IK we can estimate the magnitude of polarization of 
ferromagnets through the value of the conductance of 
the junctions. 

(2) For d-wave junctions with (110) orientation and Px- 
wave junctions, the ZBC decreases with increasing tem- 
peratures when the degree of the polarization, X, is 
small. For large X, the ZBC is an increasing function 
of T. The presence of the ZES explains these behavior. 
In d-wave junctions with (100) orientation and pj,-wave 
junctions, the ZBC is an increasing function of T inde- 
pendent of X. This is because the ZES does not appear 
at these junction interface. For px + ipy-wave junctions, 
the ZBC first decreases with increasing T then increases. 

(3) In p-wave junctions, the temperature dependence of 
the ZBC depends on the direction of the magnetization 
axis of ferromagnets because of the spin degree of freedom 
of Cooper pairs in spin-triplet superconductors. 

(4) Throughout this paper, we have calculated the ZBC 
in two ways; i) the spatial dependence of the pair po- 
tential is assumed to be the step function (non-SCF cal- 
culation), ii) the spatial depletion of the pair potentials 
are determined self-consistently (SCF-calculation). We 
have confirmed that there are no remarkable differences 
between the conductance in the non-SCF calculation and 
those in the SCF calculation. 

In this paper, effects of random potentials in ferromag- 
nets are not taken into account. Recently, there are sev- 
eral works on random scattering effects in unconventional 
superconductor junctiongSiS^iiSiiiiS. It is actually inter- 
esting to study transport properties of junctions where 
diffusive ferromagnets are attached to unconventional su- 
perconductors. In the present paper, the splitting of the 
ZBCP by magnetic fields through the Zeeman effect or 
magnetic impurities in an insulator are not taken into ac- 
coun^i&ii^i^. These are interesting and important future 
issues. 
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The present model deals with the ideal situation where 
the height of the ZBCP is only the function of the polar- 
ization of the ferromagnet. However, in real junctions, the 
interface quality and the impurity distribution near the in- 
terface may have influences on the ZBCP height. Detailed 
analysis of such factors will be discussed in coming papers. 



